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Abstract 


We calculate the heavy-light meson decay constant using lattice NRQCD 
action for the heavy quark and Wilson quark action for the light quark over a 
wide range in the heavy quark mass. Simulations are carried out on a 16^ x 32 
lattice with 120 quenched gauge configurations generated with the plaquette 
action at /3 = 5.8. For the heavy quark part of the calculation, two sets of lat¬ 
tice NRQCD action and current operator are employed. The first set includes 
terms up to Ofi/mq) both in the action and the current operator, and the 
second set up to ©(l/rug), where mq is the bare mass of the heavy quark. 
Tree-level values with tadpole improvement are employed for the coefficients 
in the expansion. We compare the results obtained from the two sets in detail 
and find that the truncation error of higher order relativistic corrections for 
the decay constant are adequately small around the mass of the b quark. We 
also calculate the IS hyperfine splitting of B meson, Mb„ — Mb and /bs/Ib 
with both sets and find that the l/rug corrections are negligible. Remaining 
systematic errors and the limitation of NRQCD theory are discussed. 
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I. INTRODUCTION 


The properties of hadrons including a heavy quark, particularly b quark, provide us 
with crucial information for constraining the Cabibbo-Kobayashi-Maskawa (CKM) mixing 
matrix of the Standard Model, which still have large uncertainties in spite of much effort 
with various approaches. For the combination the current value is 0.009 ± 0.003 

|jl[]. The large error mainly arises from uncertainties in the decay constant fs and the bag 
parameter Bb of B meson, which are needed to relate the experimentally measured — B^ 
transition rates with iV^^V^dl. It is, therefore, very important for the verification of the 
Standard Model to determine these B meson matrix elements with higher accuracy. 

The lattice technique enables us to carry out this task from the first principle of Quantum 
Chromodynamics (QCD). In this paper we concentrate on the decay constant and study 
various uncertainties in the calculation, which is also instructive for the calculation of the 
bag parameter. Extensive effort has been devoted to a lattice QCD determination of fs in 
the past [Q . We are now at the second stage where the accuracy has become the main issue. 
The largest obstacle for obtaining a reliable prediction is the large value of the b quark mass. 
A naive application of the Wilson or 0(a)-improved fermion action for the b quark causes 
a systematic error of 0{amb) in B meson quantities where mb is the b quark mass and a 
the lattice spacing. Since amb exceeds unity for lattice parameters currently accessible in 
numerical simulations, the error is expected to be large, rendering an extrapolation to the 
continuum limit unreliable. 

Non-relativistic QCD (NRQCD) is designed to remove the mass scale mq of the 
heavy quark from the theory and there are no 0{amQ) systematic errors in this approach. 
Since NRQCD is organized as a systematic f/mq expansion of the full relativistic QCD, 
relativistic errors in NRQCD are induced only by the truncation of the fjmq expansion. It is, 
therefore, possible to improve the approximation and to estimate the remaining uncertainty 
in a systematic way based on the fjmq expansion. 

Exploratory studies of the decay constant with lattice NRQCD were made by Davies et 
al. [Q and Hashimoto |^, where only a part of f/mq terms was included. A study with 
lattice NRQCD action fully including effects of the f/mq terms was carried out in Ref. 0. 
It was concluded that the magnitude of 0(l/mQ) correction is significantly larger than the 
naive expectation ~ 0{AqcD/Mp), where Mp is a pseudo-scalar meson mass, and therefore 
it was necessary to investigate the next order term in NRQCD. 

The goal of the present work is to estimate the magnitude of higher order 0(l/mQ) 
effects on the B meson decay constant. For this purpose we compare simulation results of 
the two sets of the lattice NRQCD action and the operator: the first set includes terms up 
to 0{f/mq) consistently and the second set takes into account the entire correction up to 
0{f/mq). Tree-level values with tadpole improvement are employed for the coefficients of 
the correction terms. We find that the contributions of second order in f/mq to the decay 
constants is adequately small around the b quark. In order to check the generality of the 
above statement, the IS hyperfine splitting of B meson, Mp^ — Mb and /bs/Ib are also 
investigated and similar results are obtained. Examination of systematic uncertainties other 
than the relativistic correction, such as discretization error, quenched approximation and 
the one-loop renormalization parameters, are not considered in detail in this work, leaving 
them for future studies. A preliminary report of an investigation of 0(l/mQ) corrections 
similar to our work has been reported in Ref. |^. 

This paper is organized as follows. In section we introduce the action and the current 
operators used in our calculation of the decay constants. Simulation details such as param¬ 
eter values and methods are given in Section H, followed by presentation of results for the 
decay constant and related quantities. In section 0 implications of the results and other 
systematic errors are discussed. Our conclusions are given in section ^ 
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II. LATTICE NRQCD 


NRQCD action at the tree level is obtained from a relativistic action by Foldy- 
Wouthuysen-Tani(FWT) transformation, 

C = h{il^ -mQ)h ^ ^NRQCD = (1) 

where h is a 4-component spinor of the heavy quark field and Q and y are 2-component 
fields in the NRQCD theory. The NRQCD action is represented by the following l/mq 
expansion, 


r _ p{0) I /^(l) I ^(2) I 

+ ■ ■ ■ 5 

4’= f—) (2) 

\mQj 


where the mass term is discarded since it only amounts to a constant shift in the total 
meson energy and does not affect the dynamics of the system. Lattice NRQCD action is a 
discretized version of the continuum action Wick-rotated to the Euclidean formalism. The 
discretization procedure is not unique, and we choose a form which leads to the following 
evolution equation for the heavy quark propagator: 


Gqit, x) = 0 (for f < 0) 



+^a:,o (forf>0). 


l,f) 


(3) 

(4) 


Here x = (t,x), n is the stabilizing parameter [§. Our discretization procedure is almost 
same as 


Gq{t,x) 




aHoy 
2n ) 


ul 



Gq{t-l,x), 


which was used in |0. These two discretization procedures are the best choices from the 
view of the control on the discretization error in the temporal derivative. 

Hq and 6H are defined as follows: 


H =-^ 

° 2mQ’ 

5H = J2c,5H^^\ 
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The symbols A and denote the symmetric lattice differentiation in spatial directions and 

Laplacian, respectively, and A^^^ = X]i(Af^)^. The held strengths B and E are generated 
from the standard clover-leaf operator. 

The coefficients c* in (6) shonld be determined by pertnrbatively matching the action 
to that in relativistic QCD. In the present work we adopt the tree-level valne c* = 1 for all 
i and apply the tadpole improvement to all link variables in the evolntion eqnation by 
rescaling the link variables as > U^/uq. The valne of Uq is given in Section [III A . 

The original 4-component heavy qnark spinor h is decomposed into two 2-component 
spinors Q and y after FWT transformation, 


h{x) = R 


Q(x) 



( 13 ) 


where R is an inverse FWT transformation matrix which has 4x4 spin and 3x3 color 
indices. After discretization, R at the tree level is written as follows: 


where 


R = J2R^'\ 

i 
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The tadpole improvement is also applied for these operators in onr simnlations. 

As mentioned in the Introdnction, we dehne two sets of action and FWT transformation 
{5H,R} as follows: 


where 


set I = {5Hi,Ri} and set II = {5H2,R2}, 


( 21 ) 















6Hi = 
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(22) 
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(23) 


The operators 5Hi and Ri keep only 0{l/mQ) terms while 5 H 2 and R 2 inclnde the entire 
0(l/mQ) terms and the leading relativistic correction to the dispersion relation, which is 
an 0 ( 1 / 771 ^) term. The terms improving the discretization errors appearing in Hq and time 
evolution are also included. Using these two sets, we can realize the level of accuracy of 
Oil/mq) and 0(l/mQ) for the set I and II. 


III. SIMULATIONS AND RESULTS 


A. Parameters 

Our numerical simulation is carried out with 120 quenched conhgurations on a 16^ x 32 
lattice at /? = 5.8. Each conhguration is separated by 2,000 pseudo-heat bath sweeps 
after 20, 000 sweeps for thermalization and hxed to Coulomb gauge. For the tadpole factor 
we employ uq = {PpiaqY^^ with Ppiaq the average plaquette, which takes the value uq = 
0.867994(13) measured during our conhguration generation. 

For the light quark we use the Wilson quark action with k= 0.1570, 0.1585 and 0.1600, 
imposing the periodic and Dirichlet boundary condition for spatial and temporal directions, 
respectively. The chiral limit is reached at Kc = 0.16346(7) and the inverse lattice spacing 
determined from the rho meson mass equals a~^ = 1.714(63) GeV. The hopping parameter 
Kg corresponding to the strange quark is determined in two ways from m^/rrip and rriKl'mp^ 
which yields Hg = 0.15922(39) and 0.16016(23), respectively. In our analysis we take = 
0.1600 for simplicity, except in the hnal results where the error arising from the uncertainty 

in Kg is taken into account. We use the factor y^l — ^ as the held normalization for light 
quark 0. 

For the heavy quark part of calculations, two sets of lattice NRQCD action and current 
operator are employed as described in Section O For the heavy quark mass and the stabi¬ 
lizing parameter, we use (amg, n) = (5.0,2), (2.6^, (2.1,3), (1.5,3), (1.2,3) and (0.9,4), which 
cover a mass range between 2mb and me- 

All of our errors are estimated by a single elimination jack-knife procedure. 


B. Method 

In the continuum the pseudo-scalar and vector meson decay constants are dehned by 

(0|Ao|P) = /pMp, (24) 

(OlUlU) = eJvMv, (25) 


where Aq = and U = qiih- 

The lattice counterpart is calculated in the following way. Let us dehne an interpolating 
held operator for heavy-light meson from a light anti-quark and a heavy quark held by 
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y 




( 26 ) 


Q{t,y) ^ 

V 0 y 

where Fx is the gamma matrix specifying the quantum number of the meson. The subscript 
X labels the pseudo-scalar meson (P) or the vector meson (1/), 0 is a source function and the 
superscript ‘src’ denotes the choice of smearing, i.e., ‘L’(local) or ‘S”(smeared), according 
to 

(j)^{x) = 6{x) or (jy^ix) = exp(—a|a;|^) (27) 


where a and b are hxed by a £t to the Coulomb gauge wave function measured in the 
simulation. We next dehne the local axial-vector and vector currents Jx'- 

Jx = q{x)Txh{x) 


i i 


Q{x) ^ 


(28) 


with 


rp = 7574 , ^v=li- (29) 

The inverse FWT transformation qj-q explicitly written in eq.(|^-(|T^). 

To extract the decay constant of heavy-light mesons, we calculate the following two point 
functions: 


Cy(i)=E(Ot(i.2)O?(0)>, (30) 

X 

Cn^jt) = y^((9y(t, x)C>x^(0)), (31) 

X 

c2\t)=T.{Jx{t.S)O^AQ)). (32) 

X 

We show the effective mass plot of pseudo-scalar meson at amg = 2.6 and k = 0.1585 


with local CQ^{t) and smeared source CQ^{t) for the set II in Fig.|^. We observe that the 
effective mass for the smeared source is very stable from early time slices. From inspection 
of effective mass plots such as Fig.0, we conclude that the ground state of pseudo-scalar 
meson is sufficiently isolated in the range [tminytmax] = [17)22] for both correlators, and we 
adopt this range as our fitting interval. 

Similar plots for Cjj^\t) (i=2,3,4,5) with the same parameters are shown in Fig.^ We 
find that the different operators give a consistent value for the ground state energy. Effective 
mass plots for other values of amg and k exhibit similar features. 

In order to extract the binding energy and amplitude, we fit the correlation functions to 
the following forms: 



(33) 

C'ox(^) = ^Ox^^P {-^Ox^) ’ 

(34) 

c2\t) = z2\xp(-E^j^^t). 

(36) 
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As is expected from Fig.|I|and Fig.^, all correlators with the same parameters aniQ, k and X 
give a consistent value of E irrespective of the choice of ‘L’ or ‘S” and ‘O’ or ‘ J’. In the hnal 
analysis, we £t the correlators with smeared sources to obtain Eq^ and fit other correlators 
with Eq^ as the input energy. We list results for the binding energy Eq^ in Table 

Chiral extrapolation of Eq^ is shown in Fig.^ and the extrapolated values are given in 
the last column of Table |. From Fig.^ we find that the slope become a little milder for larger 
aniQ. We also show the l/(amQ) dependence of binding energy for pseudo-scalar meson at 
K = Kc in Fig.^. 

The lattice value of the meson decay constant is obtained in terms of the amplitudes 
defined above as follows: 


_ /o yL 

i,j 

= (36) 

where and are dehned for convenience of discussions 

below. The renormalization constant Zxij representing mixing effects between different 
operators are generally needed to relate the current operator on the lattice to that in the 
continuum. A perturbative calculation for our dehnition of the NRQCD action has not yet 
been hnished, and we treat Zxij as a unit matrix in the present article. We should remark 
that the NRQCD Collaboration has recently completed such a calculation H- Their choice 
of action, however, is slightly different from ours, and the results are not applicable for our 
analysis. Thus our results are stated without the one-loop Z-factor, but incorporating the 
mean held improvement. 


C. Analysis of results 

The numerical results of for all k and Kc are tabulated in Table |1| and its 

chiral extrapolation is shown in Fig.^. We hnd from this hgure that the linear extrapolation 
is very smooth. In contrast to the binding energy, the slope tends to increase as the heavy 
quark mass becomes larger. 

Fig. H shows the l/(aMp) dependence of (/p-\/Mp)^“** at k = Kc for the set I (open 
symbols) and II (hlled symbols), where aMp is the pseudo-scalar meson mass in lattice 
units calculated as described below. Shaded bands represent the mass region corresponding 
to the B and D meson. They are estimated from the value of a~^ determined from the p 
meson mass and a string tension(~ l.SGeV) at (3 = 5.8. Solid curves are results of a fit with 
a quadratic function in l/(aMp) given by 





ai 

aMp 


+ 


02 \ 
{aMpyJ • 


(37) 


The values of the htted parameters are tabulated in Table |T|. We observe that {fp\/Mp)°° 
and ai are consistent between the set I and II within the statistical error, while 02 is different 
as expected. 

The meson mass is given by 
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(38) 


aMp ^ ZntamQ + - Eo, 

where 2 m and £”0 are the mass renormalization factor and the energy shift, respectively. 
Since perturbative results for these quantities are not fully available for our NRQCD action, 
we set aMp = amg + Eq^ in this work. For the action I, for which our one-loop results are 
available, the one-loop correction is very small (~ 5%) and we expect that effects for the 
hnal prediction for fp is not signihcant. 

It may appear at hrst sight that there are no large difference over almost all of the 
region of aMp between the results from set I (open circles) and those from set II (solid 
circles) in Fig.^ In order to investigate the differences further, we decompose the results 

into contribution of each operator 6 fp\ Fig. shows the leading contribution 6 fp'^ for 
each set. The values from the set II are larger than those from the set I, showing effects 
of the l/rriq correction in the action. In Fig.H the other current contributions 6 fp'^ (i=2 
for set I and i=2,3,4,5 for set II) are shown. As expected, the magnitude of the 0(l/mQ) 
operator(circles) is larger than the other 0(1/771^) operators(other symbols). The numerical 

data for 6 fp^ are tabulated in Table 0 

In summary, O^I/uiq) correction in the action tends to raise {fpy/MpY'^^^ while the one 
in the current lower. After all we hnd a remarkable fact that the small difference in Fig.^ 
results from a cancellation between the correction from the action and that from operators, 
and between different operators. 

In order to quantify the magnitude of 0{l/mQ) corrections, we dehne the following 
quantity. 


A/ = 


(/pVM^)f ^ - {fpVW)T 

ifpvwyr 


(39) 


where I and II corresponds to the set I and II. We show A/ in Fig.^. We see that the 
0 ( 1 / 771 ^) correction is about 3% around the B meson region, while increasing to about 15% 
around the D meson. Unless the 0(1/777^) correction in the renormalization factor is 
unexpectedly large, the smallness of 0(1/7T7 q) corrections in the region of B meson will be 
retained in the physical prediction of the decay constant. 

We have seen that (/pv^Mp)^“** does not change much with inclusion of 0(1/7?7 q) terms. 
Since this is due to a cancellation among the 0(1/7T7 q) contributions whose origin is not 
apparent, the smallness does not necessarily mean that higher order corrections of the I/ttiq 
expansion are negligible. To examine this point, we dehne the following quantity: 


|A/| ^ (I - i/SI + IVS - V<f> I + I 

+ \Sfpl\ + \Sfpl,\) IUp\fM'p)'r. (40) 


This quantity provides a conservative estimate of 0{l/mQ) correction since all 0(1/777 q) 
corrections are added. The l/{aMp) dependence of |A/| is shown in Fig.[T^ with open 
symbols, together with the result for A/ (solid symbols). If we estimate the unknown 
0(1/777 q) correction to have a magnitude ^|A/|, which is presumably an overestimation, 
we deduce that (fpy/MpYff^ would be corrected by only about 6% around B meson. On 
the other hand, there is no reason that 0(1/7?7 q) correction would be small in the D meson 
region. 

For completeness, we also show the result of the vector meson decay constant 
(/y-y/My)^®**, and the spin average and the ratio of pseudo-scalar and vector decay con¬ 
stants. 
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The numerical results for the vector meson decay constant are tabulated in Table 0. 
The results for show only small difference between the set I and II as in the 

pseudo-scalar case. Making a decomposition into current components as before, we find that 
there are cancellations among Sfy^ as in the pseudo-scalar, though to a lesser extent. The 
numerical data for the spin averaged decay constant —H^(iWW)— 

can be found in Tab.|VT[ The behavior of (/a/M)^^ as a function of is shown in 

Fig.0 where aM^y ^JaMp + 3aMy)/4:. 

Fig. shows (/p//y)^“**(circles) and (/p//y)*^^^(diamonds) with set I(open symbols) 
and set II(sohd symbols). For the numerical data, see Tab. [VII| . 


D. Other quantities 


In order to hnd how large the 0(l/mQ) corrections are in other quantities, we compare 
IS hyperhne splitting, Mp^ — Mp and fp^/fp obtained with the set I and II. Renormal¬ 
ization effects are expected to be negligible for these quantities since the light quark mass 
dependence of the renormalization constant is very small around the chiral limit. 

Fig. |T^ shows the aMp dependence of IS hyperhne splitting. The splitting linearly 
increases with l/(aMp) for both sets and 0(1/771^) terms do not affect this quantity. The 
results for aMp^ — aMp are shown in Fig.|T^. We expect from experiments that this quantity 
depends only weakly on the heavy quark mass (Mpo — Mpo = QO.lMeV and M^± — 

= 99.2MeV). Our results are consistent with this expectation including the trend that the 
mass difference increases for smaller heavy quark mass, albeit errors are large. Finally we 
show fpjfp in Fig.|^ calculated from the jackknife samples of the following ratio: 


(/pyM^)'“« \ aMp/ 


(41) 


The ratio fps/fs has phenomenological importance since it is necessary to extract the Stan¬ 
dard Model parameter \Vts\ which is up to now only poorly determined. One can see in 
Fig. [l3| - p!^ that there is no signihcant difference between the results from the two sets of 
simimtions over almost all mass region up to the D meson. Numerical results are tabulated 
in Table fTM 

With a ^ from p meson mass, we hnd Mp{amQ = 2.6) ~ 5.3GeV which is close to the 
experimental value of Mp- We therefore consider arrib = 2.6 to be the physical point for 
the b quark and convert the numerical results above from the set II into physical units. We 
obtain 


Mp* — Mp = 26 ± 9(statistical) MeV, 

Mp^ — Mp^ = 99 ± ^{statistical) ± 13{strange) MeV, 
fp 

= 1.23 ± 0.03{statistical) ± 0.03(stran(7e). 

JBd 

where ^strange’’ means the error arising from the ambiguity in Kg using m^/rrip or mp/rrip. 
The hyper hne splitting is much smaller than the experimental value of 46 MeV. It is known 
that this quantity is very sensitive to 0{a) error and quenching effects. Other quantities are 
in reasonable agreement with experiment and results of previous lattice studies. 


9 














IV. fs TO 0(l/m^) AND REMAINING SYSTEMATIC UNCERTAINTIES 


Our investigation shows that relativistic corrections of order I/toq is small in the region 

of B meson, and that higher order corrections are likely to be bound within a 5% level. 
One of the remaining source of systematic uncertainties is a discretization error of form 
0{{aKQCDY)- The leading error of this form existing in our calculation is O^aAQco) which 
appears from the Wilson fermion action since the gauge and NRQCD action have no 0{a) 
term. The characteristic size of O^aAqcD) at P=5.S is 20-30%. This error can be reduced 
to the level of 5% by the use of (9(a)-improved Wilson actions. Alternatively, one may carry 
out simulations at a larger value of (3 in order to reduce the O^aAqcD) error within the 
Wilson action for light quark. However, care must be taken in this alternative because of 
the problem of divergence of one-loop coefficient for mq < 0.6-0.8 [|^]. Such a situation 
can arise when the heavy quark mass parameter in lattice units becomes small, which will 
be encountered in simulations at large (3. These limitations in the values of (3 and l/mq 
should be kept in mind in simulations with lattice NRQCD. 

Another source of systematic errors is the deviation of the expansion parameters of 
the NRQCD action and renormalization constants of currents from their tree-level values. 
Perturbative corrections in these quantities are not negligible in general, amounting to ~ 
20% at one-loop order. The renormalization factor of the axial-vector current in the static 
limit is known to be particularly large, and l/mq corrections could also be important. We 
expect, however, that after including the one-loop correction the systematic error of this 
origin will be reduced to O(a^) ~ 5% in magnitude. 

Taking into account the uncertainties discussed above, we quote the following estimate 
from the present work for the B meson decay constant in the quenched approximation: 

fs = 184 ± 7(statistical) ± 5{relativistic) MeV, 

with discretization and perturbative errors of 20 % each. We use the static result Za = 
1 — 0.057gy{q*) |]T^ for the renormalization constant, and an average is taken of the result 
with q* = - and q* = -. 


V. CONCLUSION 

In this article we have presented a study of the 0(l/mQ) correction to the decay con¬ 
stant of heavy-light meson with lattice NRQCD and Wilson quark action in the quenched 
approximation. While the 0(l/mQ) correction to the decay constant in the static limit is 
signihcant, we hnd in our systematic study of 1 /mq expansion that the 0(l/mQ) correction 
is sufficiently small for B meson, so that there will be no need for incorporating 0(l/mQ) 
corrections unless an accuracy of better than 5% is sought for. Our examination of other 
physical quantities in the same respect also provides encouraging support to this statement. 
We have thus shown using our highly improved lattice NRQCD that the relativistic error, 
which has been one of the largest uncertainty in lattice calculations of the B meson decay 
constant, is well under control. 

Our results still have several sources of large systematic errors. In order to obtain 
with a higher precision, we need to reduce (i) the 0 (a) error by using an 0 (a) improved 
Wilson fermion action for the light quark, (ii) the 0(as) and 0(aas) errors with the fully 
one-loop corrected perturbative renormalization coefficients for both the action and the 
operator, and (iii) the scale setting and quenching error by doing simulations with full QCD 
conhgurations. The problems (i) and (ii) are currently under study, and we are planning to 
pursue (iii) soon. When these improvements are all in place, we expect to achieve a lattice 
NRQCD determination of Jb with the accuracy of less than 10%. 
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TABLES 


airiQ 

k=0.1570 

0.1585 

0.1600 

Kc=0.16346 

5.0 

0.625(7) 

0.602(8) 

0.577(12) 

0.524(16) 


0.631(7) 

0.608(8) 

0.583(11) 

0.531(15) 

2.6 

0.618(5) 

0.595(6) 

0.570(8) 

0.518(11) 


0.624(5) 

0.601(6) 

0.576(8) 

0.524(11) 

2.1 

0.613(5) 

0.590(6) 

0.565(7) 

0.512(10) 


0.618(5) 

0.594(6) 

0.569(7) 

0.516(10) 

1.5 

0.604(4) 

0.580(5) 

0.555(6) 

0.501(8) 


0.600(4) 

0.576(5) 

0.551(6) 

0.498(8) 

1.2 

0.596(4) 

0.571(5) 

0.546(6) 

0.492(8) 


0.581(4) 

0.556(5) 

0.531(6) 

0.477(7) 

0.9 

0.579(4) 

0.554(4) 

0.529(5) 

0.473(7) 


0.536(4) 

0.511(4) 

0.486(5) 

0.430(7) 


TABLE I. Binding energies of heavy-light pseudoscalar mesons for each k and set. Upper lines 
are from set I and lower lines from set II. 
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amq 

k=0.1570 

0.1585 

0.1600 

Kc=0.16346 

5.0 

0.408(13) 

0.377(13) 

0.341(14) 

0.268(17) 


0.397(12) 

0.367(12) 

0.333(13) 

0.264(16) 

2.6 

0.334(7) 

0.311(8) 

0.286(8) 

0.233(10) 


0.323(7) 

0.301(7) 

0.277(8) 

0.227(9) 

2.1 

0.310(6) 

0.289(7) 

0.266(7) 

0.219(8) 


0.298(6) 

0.278(6) 

0.257(6) 

0.212(8) 

1.5 

0.269(5) 

0.252(5) 

0.234(5) 

0.195(7) 


0.256(5) 

0.240(5) 

0.223(5) 

0.187(7) 

1.2 

0.242(4) 

0.227(4) 

0.212(5) 

0.178(6) 


0.224(4) 

0.211(4) 

0.197(5) 

0.168(6) 

0.9 

0.209(4) 

0.197(4) 

0.184(4) 

0.156(5) 


0.176(4) 

0.166(4) 

0.157(5) 

0.134(6) 

TABLE 11. Numerical results for 

and lower lines from set 11. 

for each k and set. Upper lines 

are from set I 

set 



ai 

a-2 

set I 

0.320(22) 


-0.97(11) 

0.36(10) 

set II 

0.308(20) 


-0.87(11) 

0.17(11) 


TABLE III. The coefficients obtained by fitting each data to a quadratic function. 
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amq 

Sf? 

(x 100) 

(x 100) 

5/?^ 

(x 100) 

(x 100) 

5.0 

0.290(18) 

-2.2(3) 





0.288(16) 

-2.1(3) 

-0.18(3) 

0.03(1) 

-0.066(15) 

2.6 

0.267(11) 

-3.3(3) 





0.266(10) 

-3.3(3) 

-0.46(4) 

0.13(2) 

-0.24(3) 

2.1 

0.258(9) 

-3.9(3) 





0.259(9) 

-3.9(3) 

-0.65(5) 

0.21(2) 

-0.35(3) 

1.5 

0.243(7) 

-4.9(3) 





0.252(7) 

-5.1(3) 

-1.11(7) 

0.43(4) 

-0.69(5) 

1.2 

0.235(6) 

-5.7(3) 





0.248(7) 

-6.1(4) 

-1.59(9) 

0.71(5) 

-1.07(7) 

0.9 

0.226(6) 

-6.9(4) 





0.247(7) 

-7.8(4) 

-2.78(10) 

1.35(8) 

-1.99(10) 

TABLE IV. Numerical results for 6 f^p 

at K = Kc in lattice unit, 

. Upper 

lines with set I and 


lower lines with set II. 
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amq 



(x 100) 

(x 100) 

(x 100) 

(x 100) 

5.0 

0.280(18) 

0.274(18) 

0.63(11) 





0.275(10) 

0.270(16) 

0.61(10) 

-0.16(3) 

0.003(6) 

0.017(5) 

2.6 

0.248(10) 

0.240(9) 

0.87(11) 





0.240(9) 

0.235(9) 

0.86(11) 

-0.39(5) 

0.027(9) 

0.06(1) 

2.1 

0.237(9) 

0.227(8) 

1.00(11) 





0.229(8) 

0.223(8) 

0.98(11) 

-0.53(6) 

0.05(1) 

0.09(2) 

1.5 

0.220(7) 

0.207(6) 

1.24(12) 





0.234(7) 

0.207(7) 

1.20(13) 

-0.83(8) 

0.11(2) 

0.19(2) 

1.2 

0.210(6) 

0.196(6) 

1.42(13) 





0.205(7) 

0.198(6) 

1.36(15) 

-1.11(10) 

0.19(3) 

0.30(4) 

0.9 

0.201(6) 

0.184(5) 

1.70(16) 





0.195(7) 

0.188(6) 

1.58(20) 

-1.78(16) 

0.37(5) 

0.55(7) 


TABLE V. Numerical results for {fvMy'^y'^^^ and at k = Kc in lattice unit. Upper lines 
with set I and lower lines with set II. 


16 







amg 


<5/S 

«/.S 

(x 100) 

s/S 

(x 100) 

<5/S 

(x 100) 

SIav 

(x 100) 

5.0 

0.278(17) 

0.278(18) 

-0.06(5) 





0.272(16) 

0.274(16) 

-0.07(5) 

-0.17(3) 

0.010(4) 

-0.004(2) 

2.6 

0.244(9) 

0.246(9) 

-0.18(6) 





0.237(9) 

0.242(9) 

-0.19(6) 

-0.41(4) 

0.053(5) 

-0.012(6) 

2.1 

0.232(8) 

0.235(8) 

-0.22(7) 





0.225(8) 

0.232(8) 

-0.25(7) 

-0.56(5) 

0.089(7) 

-0.018(8) 

1.5 

0.213(6) 

0.216(7) 

-0.30(8) 





0.207(7) 

0.218(7) 

-0.37(9) 

-0.90(7) 

0.193(14) 

-0.032(16) 

1.2 

0.202(6) 

0.205(6) 

-0.36(9) 





0.195(6) 

0.210(6) 

-0.49(11) 

-1.23(9) 

0.32(2) 

-0.048(25) 

0.9 

0.189(5) 

0.194(5) 

-0.46(11) 





0.180(6) 

0.203(6) 

-0.76(15) 

-2.03(15) 

0.62(4) 

-0.084(51) 


TABLE VI. Numerical results for spin averaged lattice matrix elements and their current com¬ 
ponents at K = Kc in lattice unit. Upper lines with set I and lower lines with set II. 
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amq 


S{fp/fvY^^ 

^ifp/fvY'^^ 

Sifp/fvY^^ 

Sifp/fvY^^ 


5.0 

0.954(24) 

1.055(20) 

-3.42(37) 





0.962(24) 

1.065(18) 

-3.47(35) 

1.084(58) 

14(48) 

-3.98(75) 

2.6 

0.941(24) 

1.113(16) 

-3.83(35) 





0.945(22) 

1.133(16) 

-3.89(37) 

1.186(66) 

4.8(2.1) 

-3.77(47) 

2.1 

0.925(22) 

1.136(16) 

-3.89(33) 





0.927(22) 

1.163(17) 

-4.00(37) 

1.235(67) 

4.2(1.3) 

-3.76(44) 

1.5 

0.887(21) 

1.174(15) 

-3.96(31) 





0.877(25) 

1.215(18) 

-4.24(39) 

1.335(77) 

3.92(80) 

-3.70(41) 

1.2 

0.848(22) 

1.198(16) 

-4.00(31) 





0.819(27) 

1.253(20) 

-4.45(43) 

1.426(87) 

3.80(65) 

-3.65(41) 

0.9 

0.780(24) 

1.228(17) 

-4.09(34) 





0.694(32) 

1.311(24) 

-4.92(58) 

1.56(11) 

3.61(53) 

-3.61(44) 


TABLE VII. Ratio of pseudoscalar and vector lattice matrix elements at k = Kc- Upper lines 


with set I and lower lines with set II. 
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amq 

aMp 

aMy — aMp 

(x 100) 

aMp^ — aMp 

(x 100) 

fPs/fP 

5.0 

5.524(16) 

0.81(60) 

5.31(57) 

1.271(35) 


5.531(15) 

0.86(55) 

5.27(54) 

1.258(32) 

2.6 

3.118(11) 

1.53(53) 

5.26(40) 

1.216(24) 


3.124(11) 

1.54(54) 

5.29(39) 

1.214(23) 

2.1 

2.612(10) 

1.80(53) 

5.32(37) 

1.206(21) 


2.616(10) 

1.84(55) 

5.35(36) 

1.202(21) 

1.5 

2.001(8) 

2.34(51) 

5.41(33) 

1.185(18) 


1.998(8) 

2.47(56) 

5.41(32) 

1.176(20) 

1.2 

1.692(8) 

2.81(52) 

5.48(31) 

1.171(17) 


1.677(7) 

3.05(59) 

5.47(30) 

1.159(21) 

0.9 

1.373(7) 

3.61(55) 

5.57(29) 

1.153(17) 


1.330(7) 

4.06(66) 

5.56(29) 

1.134(24) 


TABLE VIII. Simulation results at k = Kc in lattice unit. Upper lines with set I and lower 
lines with set II. 
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FIGURES 



8 12 16 20 24 

t 

FIG. 1. Effective mass plot with local source(open circles) and smeared source (solid circles) 
at mq = 2.6, K = 0.1585 with the NRQCD action including entire Xjm^Q corrections (set II). 



8 12 16 20 24 

t 

FIG. 2. Effective plot of Cj^\t) for pseudoscalar at mq = 2.6, k = 0.1585 with set II. Where 
i=2,3,4,5 correspond to circles, squares, diamonds and triangles, respectively. 
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1/k 


FIG. 3. Chiral extrapolations of of From above aniQ = 5.0(circles), 2.6(squares), 2.1 (di¬ 
amonds), 1.5(triangles), 1.2(pluses) and 0.9(crosses) in set II. 



FIG. 4. l/(amQ) dependence of chirally extrapolated binding energies of pseudoscalar(circles) 
meson from set I(open) and set Il(solid). 
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1/k 

FIG. 5. Chiral extrapolations of . From above amg = 5.0(circles), 2.6(squares), 

2.1 (diamonds), 1.5(triangles), 1.2(pluses) and 0.9(crosses) in set II. Solid lines are obtained from 
linear fits. 



FIG. 6. l/(aMp) dependence of chirally extrapolated with set I(open circles) and 

set II(solid circles). The dashed line is obtained by fitting the data from set I to quadratic function 
and the long dashed line from set II. 
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FIG. 7. l/{aMp) dependence of the leading contribution to with set I(open circles) 

and set II(solid circles). 



FIG. 8. l/(aMp) dependence of the non-leading contributions to with set I(open 

circles) and set II(solid symbols). Solid symbols 5/p^(i=2,3,4,5) are corresponding to circles, 
squares, diamonds and triangles respectively. 
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FIG. 11. l/{aMAv) dependence of spin averaged with set I(open circles) and set 

II(solid symbols). Solid line are obtained from a quadratic fit and Small symbols represent the 
extrapolated values. 



0.0 0.2 0.4 0.6 0.8 


l/(aMp) 

FIG. 12. l/{aMp) dependence of the ratio (fp/fvY^'^^ and (fp/fvY^^ at K = Kc with 
set I(open symbols) and set II(solid symbols). Circles refer to {fp/fvY°'^^i diamonds refer to 
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FIG. 13. l/{aMp) dependence of IS hyperfine splitting with set I(open symbols) and set II(solid 
symbols). 



FIG. 14. l/{aMp) dependence of Mp^ — Mp with set I(open symbols) and set II(solid symbols). 
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FIG. 15. l/{aMp) dependence of fpjfp with set I(open symbols) and set II(solid symbols). 
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